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Influence of Molecular Shape on the Dielectric 
Constant of Polar Liquids 

F. Buckley and A. A. Maryott 

The Onsager theory of the static dielectric constant of polar liquids has been modified 
to allow for the influence of molecular shape. Equations have been derived for both 
prolate and oblate spheroids in which the permanent dipole is directed along the axis of 
symmetry. These equations have been applied to a number of pure polar liquids, exclusive 
of the hydrogen bonded liquids, where reliable data were available over extended ranges 
of temperature. A satisfactory interpretation of the variation of dielectric constant with 
temperature is obtained provided the volume of the molecular cavity is maintained fixed 
and independent of temperature rather than allowed to expand in proportion to the molar 
volume as assumed by Onsager. The required molecular shapes are in general accord 
with the actual molecular geometries. The behavior of dilute solutions and mixtures with 
a nonpolar solvent are also satisfactorily correlated using the same shape parameters. It 
appears to be unnecessary to consider specific interactions between neighboring molecules 
in order to account for the static dielectric behavior of these systems. 

1. Introduction 

The Onsager [l] 1 theory has been notably successful in relating electric dipole moments 
as determined in the vapor phase with the "static" dielectric constants of polar liquids. Ex- 
ceptions are liquids in which association occurs, usually through hydrogen bonding, as in 
water and the alcohols, or for which the dipole moment is variable, as in ethylene chloride. 
For a large number of the more normal liquids, such as ethers, ketones, nitriles, and alkyl and 
aryl halides, the calculated values of dipole moment fall with more or less equal frequency 
above and below the measured gas values. The Onsager relation thus represents a rather 
satisfactory average. However, the agreement between observed and calculated values is 
far from perfect and discrepancies of the order of 10 to 20 percent are relatively common. 

Recent extensions of the dielectric theory of liquids have usually followed one of two 
patterns. In the one case, following Onsager, the environment surrounding and interacting 
with any given molecule is treated as a continuum having the macroscopic dielectric constant 
and the molecular model is altered to accommodate either optical anisotropy [2] or geometrical 
anisotropy [3, 4] or both [5, 6]. A second, more generalized and flexible approach was intro- 
duced by Kirkwood [7] and extended by others [8, 9, 10]. Interactions between the molecule 
and its immediate surroundings are treated by the specific methods of statistical mechanics 
and the medium becomes a continuum only outside of this region. However, a rigorous ap- 
plication of this procedure to liquids is usually prevented by a lack of detailed knowledge 
of the liquid structure. Consequently, it has been customary in practice to make certain 
simplifying assumptions, of a nature comparable to those inherent in the Onsager development, 
in order to obtain a "normalized" relation that should apply in the absence of specific close- 
range interactions. Experimental departures from the idealized equation are then assumed 
to be indicative of such interactions or association [2, 10, 11]. If this interpretation were 
justified, dielectric measurements would provide a very useful tool for investigating structure 
in liquids since independent evidence regarding these interactions is ordinarily lacking. 

The simplicity and explicit nature of the Onsager approach make it desirable to investigate 
its utility in more detail than has been done heretofore. In particular, it is of interest to see 
whether any reasonable and simple modification of Onsager's molecular model can account 
for the dielectric behavior of polar liquids when subjected to wide changes in temperature 
or when progressively diluted with a nonpolar solvent. Under these conditions the effect of 
molecular interactions, if important, should vary markedly. 

In this paper equations are developed for both prolate and oblate spheroids in a manner 

1 Figures in brackets Indicate the literature references at the end of this paper. 
314715—54 3 229 



formally analogous to Onsager 's treatment of the spherical case. A solution for the prolate 
spheroid comparable to that employed here has been published by Abbott and Bolton [3]. 
Analysis of a large body of data on the dielectric constants of pure liquids, mixtures, and 
dilute solutions, excluding hydrogen bonded substances, shows that deviations from the 
Onsager equation can be rather systematically interpreted in terms of the spheroidal models 
and that these shapes are at least in qualitative accord with the actual molecular geometries. 
For these systems it appears to be unnecessary, in general, to consider additional factors 
involving specific molecular interactions or dipole association. 

2. Modified Onsager Equation 

The Onsager molecular model consists of a polarizable point dipole having a moment /x 
and an isotropic polarizability a located at the center of an otherwise empty spherical cavity 
of volume v. The cavity is surrounded by a continuum having the macroscopic dielectric 
constant e. The field inside the cavity is conveniently separated into two parts, ', a reaction 

field R that is parallel to the instantaneous dipole axis and a cavity field E c parallel to the 
externally applied field E. In the absence of any applied field, the polarizable dipole is en- 
hanced by its reaction field, so that the total moment ju=/x +a#=/V[l— -(«/#) #], where 

R=gnfv and ff=(4ir/3)(2€— 2)/(2 € +l). The cavity field, E c =3e/(2e+l)E, both orients the 
permanent dipole and induces a moment that, when enhanced by its own reaction field, be- 
comes mi=aE c /[l — (a/v)g\. The average total moment in the direction of E is then m— 



m t +n cos 0, where 6 is the angle between p and E. This model leads to the following relation 
between the dielectric constant and the molecular quantities : 

(€-l)(2e-l) v _4wN r a 
9e ° 



[ <* | Mo 

i-S. •w(i-5»)'_ 



(i) 



where k is Boltzmann's constant, T the absolute temperature, V the molar volume, and N 
Avogadro's number. The left-hand side of this expression is often called the Onsager-Kirkwood 
polarization, P. 

It is convenient to express differences between theory and experiment in terms of a 
deviation factor, G, defined as Mo(caic.)/Mo(gas). A value of unity for G indicates agreement 
between the dipole moment calculated from the dielectric constant of the liquid and the 
measured gas value. The values of G will depend upon the choice of polarizability and the 
volume assumed for the cavity. Onsager chose the cavity volume to be the same as the average 
volume occupied by the molecule, i. e., v=V/N, since eq (1) then reduces to the Clausius- 
Mossotti equation for a nonpolar liquid, a relation which, though not exact, is a good approxi- 
mation in such cases. Although the polarizability should include a small "atomic" contribu- 
tion as well as the optical part, the former is seldom known with any degree of accuracy for 
polar substances even from measurements on the vapor and will be neglected in the analyses 
to follow. Values of G derived from eq (1), hereafter designated as GWager (#on.), have been 
obtained by using the relations 4irNa/3 = R D (molar refraction for the Na D line), and 

, 3R D 3(^-1) 

a rv= — = — ^— = — • 

1 4tV 4x«+2) 

If one alters the Onsager model only to the extent of replacing the spherical cavity by one 
of equal volume having the shape of either a prolate or an oblate spheroid with the permanent 
dipole directed along the axis of symmetry, one obtains an equation of the form 



(e-l)(26+l) „_47riV 
9e 3 



[ oka I Mo 
l-pk B SkT(l~gk B y_ 



(2) 
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Figure 1. Reaction field correction coefficient as a Figure 2. Cavity field correction coefficient as a 

function of dielectric constant for various values of function of dielectric constant for various values of 

molecular eccentricity. molecular eccentricity. 

The coefficients k a , k Rl and k c depend on the dielectric constant of the medium and the shape 
of the spheroid. k R is the ratio of the reaction field of the spheroid to that of the equivalent 
sphere, and k c is the corresponding ratio for the cavity fields. For the sphere, k a =k R =k c =l, 
and eq (2) reduces to Onsager's equation. Details of the mathematical development, as well 
as the explicit functions defining k R , k c , and k a , are found in the appendix. To facilitate 
application of eq (2), it has been assumed that the product k a k c =l. This approximation is 
adequate for the present purpose because k a k c is generally close to unity and practically inde- 
pendent of temperature. 

Figure 1 shows the variation in reaction-field coefficient k R with dielectric constant for 
both prolate and oblate spheroids of various eccentricities e. A similar plot of the cavity- 
field coefficient k c is given in figure 2. For prolate spheroids, where k R and k c are both less 
than unity, the Onsager relation should give G factors less than unity. In like manner, oblate 
spheroids should lead to values of Gon. greater than unity because both k R and k c are greater 
than unity. 

3. Analysis of Data on Pure Liquids 

Initial attempts to interpret deviations from Onsager's equation in terms of eq (2) for a 
number of pure polar liquids over rather extended ranges of temperature were only partially 
successful. In particular, G factors adjusted to unity at some temperature by suitable choice 
of shape and eccentricity did not remain substantially independent of temperature. These 
variations were not systematically correlated with any simple picture of dipole association. 
The difficulty appeared to be inherent in Onsager's assumption requiring the cavity volume 
to vary with the thermal expansion of the liquid and was largely overcome when the cavity 
volume was allowed to remain independent of temperature. Moreover, analyses made by the 
authors on a number of nonpolar liquids, as well as the published work of others [12], show 
that the fit of eq (1) is improved if this volume is maintained constant. Although the empirical 
evidence for a constant cavity volume is sufficiently convincing, it should be emphasized that 
in the electrostatic theory this equivalent volume is determined by the electrostatic properties 
of the molecule and might well be essentially independent of the temperature. 
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G factors for a number of pure polar liquids over a wide range of temperature, together 
with other pertinent data, are given in table 1. 



Table 1. Summary of the analysis of the dielectric data for pure liquids' 1 



Substance 



T(°K) 



RdIV 



G0n. 



Gon. 



Substance 



T(°K) 



Fd/V 



Go n . G*On. 



PROLATE 



CH3CL.-. 
mo=1.87»>_ 

€=0.50.- 

CH 3 Br... 
/uo=1.80_. 
Rd = 14.6 
€ = 0.71.. 

CH3I 

/xo=1.65-. 
i?z> = 19.3 
€=0.85.- 



C 2 H 5 Br 

mo=2.03,-- 
Rd = 19.1 - 
€ = 0.55---- 



C2H5I 

M0=L91-. 
Bd = 24.3 
6=0.72.- 

CaH 5 CN_. 
/uo=4.03-. 
Rd = 15.8 
€ = 0.73.-. 

W-C4H9CL 
M0=2.05-. 
Rd = 25 A 
e=0.40_- 



203 


16.9 


0. 251 


87 


0.89 


1.00 


213 


15.8 


.248 


88 


.89 


1.00 


223 


14.9 


.244 


88 


.89 


0.99 


233 


14.0 


.240* 


88 


.88 


.99 


243 


13.3 


.236 


89 


.88 


.99 


253 


12.6 


.233 


90 


.88 


.99 


173 


16.9 


.305 


70 


.75 


1.01 


193 


14.9 


.297 


72 


.75 


1.00 


213 


13.3 


.289 


73 


.75 


1.00 


233 


11.9 


.281* 


75 


.75 


.99 


253 


10.8 


.274 


76 


.75 


.99 


273 


9.82 


.266 


78 


.75 


.99 


213 


9.74 


.340 


57 


.59 


1.00 


233 


8.87 


.332* 


59 


.59 


1.00 


253 


8.13 


.325 


60 


.59 


.99 


273 


7.53 


.317 


62 


.60 


1.00 


293 


7.00 


.309 


64 


.60 


1.00 


313 


6.50 


.301 


65 


.60 


.99 


173 


17.1 


.297 


77 


.86 


.99 


193 


15.2 


.290 


79 


.87 


1.00 


213 


13.6 


.283 


81 


.87 


1.00 


233 


12.3 


.277 


83 


.87 


1.01 


253 


11.2 


.270 


84 


.88 


1.01 


273 


10.2 


.263 


86 


.88 


1.01 


293 


9.31 


. 256* 


87 


.87 


.99 


313 


8.44 


.248 


87 


.86 


.98 


193 


11.6 


.336 


62 


.69 


.95 


213 


10.6 


.329 


65 


.70 


.96 


233 


9.74 


.322 


66 


.71 


.96 


253 


9.01 


.315 


69 


.72 


.98 


273 


8.38 


.308 


71 


.73 


.99 


293 


7.82 


.300* 


74 


.74 


1.00 


313 


7.29 


.293 


75 


.74 


1.00 


333 


6.80 


.286 


77 


.75 


1.00 


273 


31.0 


.230 


75 


.76 


1.01 


293 


27.2 


.224* 


73 


.73 


.97 


323 


24.3 


.216 


76 


. 75 


.99 


183 


12.2 


.276 


84 


.92 


1.00 


203 


11.0 


.269 


87 


.93 


1.01 


223 


9.94 


.263 


88 


.93 


1.01 


243 


9.07 


.257 


90 


.93 


1.01 


263 


8.32 


.251 


92 


.94 


1.01 


283 


7.66 


.246 


93 


.94 


1.01 


303 


7.09 


.240* 


94 


.94 


1.01 


323 


6.56 


.234 


95 


.94 


1.01 


343 


6.05 


.227 


96 


.93 


1.00 



CH3N02— 

M0=3.44_. 
R D = 12.5 
€ = 0.31... 

CH2CI2--- 
juo=1.56- . 
7?z> = 16.6 
€ = 0.60.-. 



CHCI3 — -. 
/*o=1.01... 
Rd = 21.4 

6 = 0.81.-. 



«-C 4 H9Br. 

Ain=2.08-. 
/?d = 28.3 
6=0.65.- . 



W-C4H9I--. 

M0=2.12_. 
/?o=33.5 

6=0.87--. 



CeHsCL-- 

mo=1.70-. 

/?D=31.1 

6=0.74.. 

CeHsBr-.. 

M0 = 1.70_. 
#D = 34.0 
6=0.82-_. 

C6H5NH2- 

mo=1.53.. 
#0=30.6 
6=0.27.. 



C 6 H 5 N0 2 - 
M0 = 4.25_. 
Rd = 32.7 
6=0.62.-. 



C 6 H 5 CN-.. 
M o = 4.42_. 
7?d = 31.6 
6—0.76... 



OBLATE 



183 


11.1 


.296 


.71 


.78 


203 


10.1 


.290 


.74 


.79 


223 


9.26 


.284 


.76 


.80 


243 


8.52 


.278 


.78 


.81 


263 


7.88 


.272 


.80 


.82 


283 


7.32 


.266 


.82 


.83 


303 


6.80 


.260* 


.83 


.83 


323 


6.35 


.254 


.86 


.83 


343 


5.93 


.248 


.86 


.83 


363 


5.54 


.242 


.87 


.83 


193 


8.89 


.325 


.54 


.59 


213 


8.18 


.319 


.56 


.60 


233 


7.53 


.313 


.57 


.60 


253 


7.00 


.306 


.59 


.61 


273 


6.54 


.300 


.61 


.62 


293 


6.12 


.294* 


.62 


.62 


313 


5.74 


.288 


.63 


.62 


333 


5.42 


.282 


.65 


.63 


353 


5.11 


.276 


.66 


.63 


223 


7.28 


.327 


.68 


.72 


243 


6.76 


.321 


.70 


.73 


263 


6.26 


.315 


.71 


.73 


283 


5.82 


.309 


.72 


.73 


303 


5.46 


. 303* 


.73 


.73 


323 


5.15 


.297 


.75 


.73 


343 


4.89 


.291 


.76 


.74 


373 


4.52 


.282 


.78 


.74 


403 


4.17 


.272 


.79 


.73 


274 


5.74 


.330 


.64 


.65 


298 


5.39 


.322* 


.66 


.66 


313 


5.18 


.318 


.67 


.66 


328 


4.96 


.314 


.68 


.66 


273 


7.38 


.341 


.94 


.95 


283 


7.13 


.339 


.94 


.95 


293 


6.89 


.336* 


.95 


.95 


303 


6.66 


.333 


.96 


.95 


313 


6.43 


.330 


.96 


.94 


323 


6.22 


.327 


.96 


.94 


283 


37.6 


.322 


.97 


.98 


303 


33.9 


.317* 


.96 


.96 


323 


30.6 


.312 


.95 


.93 


343 


27.6 


.306 


.94 


.90 


363 


24.9 


.301 


.92 


.87 


383 


22.7 


.296 


.91 


.85 


403 


20.8 


.290 


.91 


.84 


423 


19.3 


.285 


.91 


.83 


473 


16.0 


270 


.91 


.80 


273 


27.6 


.313 


.64 


.65 


291 


25.9 


.309* 


.65 


.65 


313 


24.0 


.303 


.68 


.67 


323 


23.3 


.300 


.69 


.67 


333 


22.7 


.297 


.70 


.68 


343 


22.1 


.295 


.72 


.69 



283 


39:1 


0.236 


1.03 


1.05 


1.00 


303 


35.9 


.230* 


1.05 


1.05 


1.00 


323 


32.9 


.224 


1.07 


1.05 


1.00 


343 


30.1 


.218 


1.08 


1.05 


1.00 


363 


27.6 


.213 


1.09 


1.03 


1.00 


173 


17.0 


.301 


1.09 


1.23 


1.00 


193 


15.0 


.294 


1.11 


1.22 


1.00 


213 


13.4 


.287 


1.13 


1.22 


1.00 


233 


12.0 


.280 


1.15 


1.22 


1.01 


253 


10.9 


.273 


1.17 


1.22 


1.01 


273 


9.92 


.266 


1.19 


1.22 


1.01 


293 


9.09 


.259* 


1.21 


1.22 


1.01 


313 


8.27 


.252 


1.21 


1.20 


1.00 


213 


6.76 


.294 


1.46 


1.56 


1.00 


233 


,6.12 


.287 


1.46 


1.53 


1.00 


253 


5.61 


.280 


1.47 


1.52 


1.00 


273 


5.22 


.274 


1.50 


1.52 


1.02 


293 


4.81 


.267* 


1.49 


1.49 


1.01 


313 


4.47 


.260 


1.50 


1.48 


1.01 


333 


4.12 


.253 


1.47 


1.43 


1.00 



(CH 3 )3CC1 
mo = 2.13_- 
7?z> = 25.9. 
6 = 0.67.. . 

(CH 3 )20— . 
mo = 1.30-. 
i?z) = 13.3- 

6 = 0.58..- 



(C 2 H 5 )20... 
mo = 1.15-. 
Rd = 22. 5 
6=0.89... 



.97 
.99 
1.00 
1.00 
1.01 
1.01 
1.01 
1.01 
1.01 
1.01 



.99 
1.00 
1.00 
1.00 
1.00 
1.00 
1.00 
1.00 

1.01 
1.01 
1.01 
1.00 



1.00 
1.01 
1.01 
1.01 

1.00 
1.00 
1.00 
1.00 
.99 



1.25 
1.23 
1.19 
1.16 
1. 12 
1.09 
1.06 
1.05 
1.01 

0.97 



1.00 
1.01 
1.03 



253 


12.2 


.250 


1.31 


1.36 


263 


11.5 


.246 


1.31 


1.35 


283 


10.4 


.239 


1.32 


1.34 


303 


9.37 


.233* 


1.32 


1.31 


298 


5.02 


.189* 


1.16 


1.16 


318 


4.54 


.179 


1.17 


1.15 


338 


4.07 


.169 


1.17 


1.13 


358 


3.59 


.156 


1.17 


1.10 


378 


3.12 


.137 


1.19 


1.10 


398 


2.37 


.102 


1.20 


1.08 


233 


5.86 


.237 


1.63 


1.71 


253 


5.30 


.230 


1.62 


1.67 


273 


4.80 


.224 


1.60 


1.63 


293 


4.34 


.217* 


1.57 


1.57 


313 


3.97 


.209 


1.55 


1.53 


333 


3.67 


.202 


1.55 


1.51 


353 


3.40 


.194 


1.52 


1.48 


373 


3.14 


.185 


1.53 


1.46 



1.01 
1.01 
1.00 
0.99 

1.03 
1.03 
1.01 
0.99 
1.01 
1.00 

1.00 
1.01 
1.01 
0.99 
.99 
1.00 
1.00 
1.01 
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Table 1. Summary of the analysis of the dielectric data for pure liquids — Continued 



Substance 



T(°K) e F Dl 'V Gon. Chn. 0* 



Substance 



T(°K) t Fd/V C on. Gon. 0« 



O B L AT E— Continued. 



(CH 3 ) 2 CO 

M0=2.88 

7?d = 16.2 

e=0.55 

C2H5COCH3. 

M o=(2.88)d ._. 

F D = 20.7 

e=0.55 

(C 2 H 5 )2CO«.... 
Mo=(2.88)d— . 

7?z> = 25.3 

6=0.60 

rc-C 3 H7COCII 3 
M0=(2.88)d._. 

Rd = 25.4 

6=0.50 



193 


34.5 


.251 


1.03 


1.11 


0.96 


213 


31.3 


.245 


1.07 


1.14 


.98 


233 


28.4 


.239 


1.10 


1.15 


1.00 


253 


25.9 


.233 


1.13 


1.17 


1.02 


273 


23.7 


.227 


1.16 


1.18 


1.02 


293 


21.5 


.220* 


1.18 


1.18 


1.02 


313 


19.4 


.214 


1.18 


1.17 


1.01 


213 


27.2 


.255 


1.10 


1.17 


1.00 


233 


24.6 


.249 


1.13 


1.18 


1.01 


253 


22.3 


.243 


1.15 


1.19 


1.02 


273 


20.3 


.237 


1.18 


1.19 


1.02 


293 


18.5 


.231* 


1.19 


1.19 


1.02 


303 


16.8 


.225 


1.20 


1.18 


1.02 


333 


15.3 


.219 


1.20 


1.17 


1.01 


353 


13.9 


.213 


1.20 


1.15 


1.00 


273 


18.9 


.245 


1.26 


1.28 


1.05 


293 


17.0 


.239* 


1.26 


1.26 


1.04 


313 


15.3 


.233 


1.25 


1.23 


1.02 


333 


13.8 


.227 


1.24 


1.21 


1.00 


353 


12.5 


.221 


1.23 


1.18 


0.98 


373 


11.5 


.215 


1.24 


1.17 


.98 


213 


22.0 


. 201 


1.04 


1.10 


.97 


233 


20.2 


.255 


1.08 


1.13 


1.00 


253 


18.4 


.249 


1.10 


1.13 


1.01 


273 


16. , r > 


.243 


1.13 


1.14 


1.01 


293 


15.5 


.238* 


1.15 


1.15 


1.02 


313 


14. 1 


.232 


1.16 


1.14. 


1.01 


333 


12.9 


.226 


1.16 


1.13 


1.01 


353 


11.7 


.220 


1.16 


1.12 


0.99 


373 


10.8 


.212 


1.19 


1.12 


1.00 



i-CMIoCOCHs- 

mo = (2.88) d ._. 

F D = 30.3 

6 = 0.40 

(n-C 3 H 7 )2CO... 

mo = (2.88) d ._.. 

/?/) = 34.4 

6=0.56 

B-CgHnCOCHj 

M0=(2.88)d... 

/?p = 35.2 

6 = 0.40 



n-CeHisCOClh 

M0=(2.88)d __ 

/?/) = 40.3 

6 = "---- 
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17.4 


.260 


1.06 


1.11 


253 


15.9 


.254 


1.09 


1.12 


273 


14.5 


.249 


1.11 


1.12 


293 


13.1 


.243* 


1.11 


1.11 


313 


11.8 


.237 


1.09 


1.08 


333 


10.7 


.231 


1.09 


1.06 


353 


9.75 


.225 


1.08 


1.04 


373 


8.90 


.220 


1.08 


1 . 02 


253 


15.1 


.256 


1.15 


1.18 


273 


13.8 


.251 


1.17 


1.18 


293 


12.6 


.246* 


1.18 


1.18 


313 


11.4 


.240 


1.17 


1.15 


333 


10.4 


.235 


1.16 


1.13 


353 


9.46 


.230 


1.17 


1. 11 


373 


8.65 


.223 


1.15 


1.09 


393 


8.00 


.217 


1.15 


1.08 


253 


14.3 


.263 


1.06 


1.09 


273 


13.1 


.257 


1.08 


1. 10 


293 


12.0 


.252* 


1.09 


1.09 


313 


10.9 


.246 


1.08 


1.07 


333 


9.93 


.241 


1.08 


1.05 


353 


9.08 


.235 


1.08 


1.04 


373 


8.27 


.229 


1.07 


1.01 


393 


7. 61 


.223 


1.06 


1.00 


413 


7. 10 


.217 


1.08 


1.00 


253 


12.5 


.269 


1.01 


1.04 


273 


11.5 


.263 


1.02 


1.04 


293 


10.4 


.258* 


1.02 


1.02 


313 


9.42 


.252 


1.01 


1.00 


333 


8.70 


.247 


1.02 


0.99 


353 


8.01 


.241 


1.02 


.98 


373 


7.42 


. 236 


1.02 


.97 


393 


6.90 


.230 


1.03 


.96 


413 


6.49 


.224 


1.04 


.97 


433 


6. 10 


.218 


1.06 


.97 



1.03 
1.04 
1.C4 
1.03 
1.00 
0.98 
.97 
.95 

1.01 
1.01 
1.01 
1.00 
0.99 
.97 
.95 



1.01 
1.02 
1.01 
0.99 
.98 
.97 
.95 



1.04 
1.04 
1.02 
1.00 
0.99 
.98 
.97 
.96 
.97 
.97 



a Sources of data: A. A. M.uyoti and E. R. Smith, Table of dielectric constants of pure liquids, NBS Circular 514 (1951); A. A. Maryott, 
and Floyd Buckley, Table of dielectric constants and electric dipole moments of substances in the gaseous state, NBS Circular 537 (1953). 

b All values of dipole moment given in Debye units (1 Debye = 10 -18 esu). 

°Data at the lower temperatures have been omitted since the tabulated values of dielectric constant are inconsistent with the graphical repre- 
sentation (R. N. Cole, .1. Cnem. Phys. 9, 251 (1961)). 

d Value assumed to be the same as for acetone. 



The factors designated as G% n . and G e were obtained from eq(l)and (2), respectively, with the 
assumption of a fixed cavity volume. In these cases, the value of B D jV noted by an asterisk, 
a value usually pertaining to room temperature or its vicinity, was rather arbitrarily used. 
Estimates of the parameter a/v are not critical, and values that differ reasonably from the 
selected values do not entail significant changes in the results of the calculations or the inter- 
pretation of the data. Compounds have been classified as prolate or oblate, depending upon 
whether 6?g n . is less or greater than unity. Using the reaction and cavity field coefficients 
appropriate to this shape, values of eccentricity were found that fitted eq (2). As the cor- 
responding deviation factors G e have been adjusted near unity, they are only indicative with 
respect to their constancy with temperature and with respect to any correlation between the 
required cavity shape and the molecular geometries of the simpler types of molecules considered. 
Table 1 shows that the values of G 0n . frequently increase with rising temperature. Such 
behavior, if confined to cases where 6? 0n . is less than unity, might be rather simply explained 
in terms of molecular interactions. Thus for molecules simulating prolate spheroids, dipole 
interactions would be such as to favor antiparallel alinement [13] and thereby tend to reduce 
the dipolar polarization. Value of G 0n . should then be less than unity but approach unity at 
more elevated temperatures as a consequence of the diminishing effect of such interactions. 
I n like manner, values of G 0n . greater than unity, if interpretable in terms of parallel association 
in conformity with the general oblate geometry of this group, should also approach unity with 
rising temperature. However, except for ethyl ether, no such trend is evident. Furthermore, 
the values of G* n., calculated on the basis of a fixed cavity volume, generally show little or no 
variation with temperature for the prolate cases. Thus the evidence for association in these 
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substances would appear to depend upon Onsager's assumption of an expanding cavity volume. 

With the exception of nitrobenzene and several of the ketones, the modifications for cavity 
shape and cavity volume lead to good agreement with the experimental data. The values of 
G e usually vary from unity by no more than 1 or 2 percent over wide ranges of temperature, 
180° in some instances, and are probably commensurate with the experimental uncertainties 
in the values of dielectric constant. The series CH 2 C1 2 , CHC1 3 , and (C 2 H 5 ) 2 is particularly 
indicative. Although the values of G 0n . are all greater than unity, they show a different type 
of dependence on temperature in each case. These differences are readily interpre table in 
terms of the modified relation. For CH 2 C1 2 , which is typical of substances having relatively 
high values of dielectric constant and only moderate values of eccentricity, G 0n . shows a pro- 
nounced increase with rising temperature. It is evident from figures 1 and 2 that in this case 
the correction factors for the reaction and cavity fields, k R and k C) are insensitive to the changes 
in dielectric constant accompanying changes in temperature. As G e is independent of tem- 
perature, the variation in 6r n- arises almost entirely from the assumption that the cavity 
volume expands with the thermal expansion of the liquids. The values of G 0n . are relatively 
independent of temperature for CHC1 3 and decrease with increasing temperature for (C 2 H 5 ) 2 0. 
As these substances have unusually low values of dielectric constant and require exceptionally 
high values of eccentricity, k R and k t now vary significantly with temperature. 

Nitrobenzene fits the Onsager relation more closely than most substances, although there 
is a definite decrease in G 0n . with increasing temperature. The modified relation leads to a 
greater dependence of G on the temperature. On the presumption that this variation is 
attributable to association, the shape was chosen to fit the data at the highest temperature. 
The value of eccentricity required on this basis is in reasonable accord with that anticipated 
in comparison with chlorobenzene and other benzene derivatives of similar shape. 

Of the eight ketones listed only acetone, methyl ethyl ketone, and methly 7i-propyl ketone 
conform reasonably well to eq (1). Diethyl ketone and higher members of the series show a 
definite downward trend in G e with increasing temperature. This trend can hardly be attributed 
to association in view of the normal behavior of the lower members for which such association 
should be most pronounced. However, all of the ketones except acetone have flexible alkyl 
groups and it is not unreasonable to suppose that the effective shapes, especially of the higher 
members, may vary with the temperature. The observed trend in G e then suggests a change 
to a less oblate or more prolate shape as the temperature increases. Although similar com- 
plications might conceivably arise because of sufficient molecular flexibility in several other 
cases, namely, ethyl ether and the ^-butyl halides, only the higher ketones appear to exhibit 
any abnormal behavior. 

Table 2 contains a number of additional liquids for which reliable data were not available 
over a sufficient temperature range to be included in table 1. The values of eccentricity re- 
quired to fit eq (2) in the vicinity of room temperature are listed in the final column. 

It is of interest to see how the shapes derived from the dielectric analyses correspond to es- 
timates of the molecular shapes for those compounds where the dipole moment is directed along 
an axis of symmetry. In table 3 values of the axial ratio bja corresponding to the listed ec- 
centricities are compared with those obtained by scaling models constructed from the Fisher- 
Hirshfelder-Taylor atomic models. The dimension a refers to the dipole axis while b was taken 
as the largest dimension perpendicular to this axis. Some molecules containing flexible ethyl 
groups are included. For these cases, the values of b were obtained from the most extended 
configuration and, consequently, these axial ratios represent upper limits. In no case do the 
two procedures lead to a definite conflict in the assignment of oblate and prolate character. 
However, the molecular geometry of t-bntyl chloride is virtually spherical, whereas eq (2) 
requires a decidely oblate shape. A similar discrepancy, but in the reverse sense, is noted for 
methylene bromide. On the other hand, the high eccentricities needed for ethyl ether and tri- 
ethylamine are close to those estimated for the most extended configurations. In the prolate 
cases the agreement is consistently good, except for iodobenzene. 
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Table 2. 



Summary of the analysis of the dielectric data for pure liquids for which data over an extended temperature 

interval are not available a 



Substance 



A 


Rd 


T 


e 


Rd/V 


Gon. 



PROLATE 



CH 3 CN.... 
w-CsHtCN. 
W-C4H9CN. 
C2H5NO2-- 
W-C3H7NO2 

71-C3H7I--. 

Ti-CaHyBr.. 
n-CsHuBr.. 
w-C 6 Hi 3 Br_. 
7i-C 7 Hi5Br.. 

n-CsHnBr.. 
TO-C 9 Hi 9 Br_. 
w-CioH 2 iBr. 
n-Ci2ll25Br. 
w-CuH^Br. 

C 6 H 5 F 

C 6 H 5 I 

CH 2 Br 2 



3.97 


11.1 


293. 


37.5 


0.212 


0.80 


4.07 


20.4 


293 


20.3 


.234 


.60 


4.12 


25.0 


293 


17.4 


.239 


.58 


3.62 


17.0 


303 


28.1 


.235 


.96 


3.66 


21.4 


303 


23.2 


.239 


.95 


2.04 


28.9 


293 


7.00 


.297 


.68 


(2. 10) 


26.6 


298 


8.09 


.259 


.84 


(2. 10) 


33.1 


298 


6.31 


.265 


.82 


(2. 10) 


37.7 


298 


5.82 


.266 


.83 


(2. 10) 


42.1 


298 


5.33 


.267 


.81 


(2. 10) 


46.7 


298 


5.00 


.268 


.81 


*(2. 10) 


51.5 


298 


4.74 


.270 


.81 


(2. 10) 


56.4 


298 


4.44 


.271 


.79 


(2. 10) 


66.1 


298 


4.07 


.272 


.78 


(2.10) 


74.6 


298 


3.84 


.272 


.78 


1.59 


26.0 


293 


5.47 


.277 


.83 


1.70 


39.3 


293 


4.63 


.353 


.45 


1.4J 


21.9 


283 


7.77 


.320 


.99 



OBLATE 



CH3CHCI2— 
CH3CCI3-— . 

(CH 3 )2CC1 2 -. 
CH 3 COC 6 H 5 . 
(n-C 3 H 7 )20--. 

(n-C4H 9 )20--. 

(C 2 H 5 ) 2 S 

(C 2 H 5 )3N_... 



2.06 


21.2 


291 


10.0 


0.252 


1.03 


1.79 


26.2 


293 


7.52 


.263 


1.12 


2.27 


25.8 


293 


10.2 


.248 


1.09 


3.02 


34.5 


298 


17.4 


.294 


1.12 


1.21 


31.7 


299 


3.39 


.230 


1.16 


1.17 


40.8 


298 


3.06 


.239 


1.15 


1.54 


28.5 


298 


5.72 


.263 


1.17 


0.66 


33.8 


298 


2.42 


.242 


1.31 



0.67 
.90 
.91s 
.30 
.30 

.80 
.60 
.65 
.65 
.67 



.68 
.71 
.72 
.72 

.64 
.95 



.50 
.45 
.45 
.65 

.62 
.58 



a Sources of data: A. A. Maryott and E. R. Smith, Table of dielectric constants of pure liquids, NBS Circular 514 (1951); A. A. Maryott and 
Floyd Buckley, Table of dielectric constants and electric dij>ole moments of substances in the gaseous state, NBS Circular 537 (1953). 
b All values of dipole moment given in Debye units (1 Debye=10- 18 esu). 



Table 3. — Comparison of axial ratios determined from dielectric data and from molecular geometry 



Substance 


bja (dielectric) 


b/a (geometry) 


Substance 


b/a (dielectric) 


b/a (geometry) 


PROLATE 


OBLATE 


CH3CI 


0.9 
.7 
.5 
.7 
.6 
.7 
.6 
.3 
.7 


0.8 


CH2CI2 


1.3 
1.0 
1.7 
1.2 
1.1 
1.4 
1.2 
2.2 
1.2 
1.4 
1.2 
1.7 
1.0 


1.4 
1.7 
1.4 
1.0 
1.2 
1.0 
1.7 
2.0 
1.3 
1.8 
1.8 
1.8 
0.9 


CH 3 Br 




7 
7 
7 
8 
8 
8 
7 
7 


CH 2 Br 2 

CHCI3 

CH3CCI3 


CH3L 


CH3CN 


C 6 H 5 N0 2 


(CH 3 ) 2 CC1 2 

(CH 3 ) 3 CC1 

(CH 3 ) 2 


C 6 H 5 C1 


C 6 H 5 Br 


C 6 H 5 I 

C 6 H 5 CN 


(C 2 H 5 )20 


(CH 3 ) 2 CO 






(C 2 H 5 ) 2 CO 

(C 2 H 5 )2S... 


(C 2 H 5 ) 3 N 






CH3NO2 



Variations in the values of eccentricity obtained from the dielectric analysis are usually in 
the direction anticipated for members of homologous series and other related compounds. For 
example, the prolate eccentricity increases with increasing size of the halogen atom, 
F<Cl<Br<I, for the alkyl and aryl halides. The prolate eccentricity also increases with 
increasing chain length in homologous series of alkyl halides and nitriles, excepting certain 
first members of the series. 
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4. Analyses of Data on Solutions and Mixtures 

On the assumption that the Onsager-Kirkwood polarizations conform to the additive 
mixture law, eq (1) becomes for a two component mixture 



~12 — 7, -^12=71^1-1-/2^2 K~ 71 

\f€\ 2 6 



a 



_i-a* a*r(i-i».) 



4ttA^ 



/* 



M2 



_!-**, 3.T(l-^ 12 ) 2 



(3) 



where / is the mole fraction, and P x and P 2 are the polarizations of the components in the 
mixture. The subscripts 1, 2, and 12 refer to the two components and the mixture, respectively. 
For dilute solutions of a polar solute in a nonpolar solvent, it is convenient to expand eq (3) 
in the Hedestrand [14] manner. The polarization of the solute at infinite dilution P^° is then 
given by 



p;=p 1 g+p I 



2e?+l 



8x- 



Vi 



e,(2«, + !)(«,- 1) 



Va/ 2 A-o d t \A/,A-o' 



(1-1^(2.4-1)^ 

where M is the molecular weight, d the density, and 

p _ (ei-l)(2ei+l) Mi 
1 9ei d 1 ' 

If the Clausius-M ossotti relation is accepted for the pure nonpolar solvent so that 

3(6l~l) 



«iM= 



4^(61 + 2)' 



the second term in brackets simplifies to 



Vi 



6Pid, 



(l~^i)(2ei + ir 



Mi(2ei + 1) 2 



For spheroidal molecular shapes, 
™ 47riV 



V K 



1 ^2 j 



3kT 



0-2*-) 



(5) 



where the coefficients k R and k c are functions of the eccentricity of the solute and the dielectric 
constant of the solvent. 

Data on six polar solutes in a number of nonpolar solvents having dielectric constants in 
the range 1.9 to 2.8 are given in table 4. As in the case of the pur 3 liquids, values of 6? 0n .fall 
below unity for C 2 H 5 Br, C 6 H 5 Br, CH 3 CN, and C 6 H 5 CN, and above unity for (CH 3 ) 3 CC1 and 
CHC1 3 , indicating that these differences may be again correlated with molecular shape. When 
the data are interpreted in terms of the modified Onsager expression (eq 5), using, values of 
k R and k c corresponding to the eccentricities derived from the pure liquids, the deviation 
factors G e are usually much closer to unity. The correction for shape usually leads to a value 
of the dipole moment within several percent of the measured gas value. Exceptions are 
CH 3 CN and CHCI3, where the corrections are decidedly too small. 
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Figure 3. G- factor versus mole-fraction plots of mixture at 20° C. 

X, GonsaEer, Of Q»'rftt mole fraction of polar constituent. 

a, A. Audsley and F. R. Goss, J. Chem. Soc. 1941,864;,. b, A. Audsley and F. R. Goss, J. Chem. Soc. 1942,358: c, A. Audsley and F. R. Goss, 
J. Chera. Soc. 1942, 497; d, F. R. Goss, J. Chem. Soc. 1937, 1918; e, D. P. Earp and S. Glasstone, J. Chem. Soc. 1935, 1709;f, G. Brtegleb, Z, 
physik.Chem. B16, 249 (1932); g, D. D. Hammick, A. Norris, and L. E. Sutton, J. Chem. Soc. 1938, 1755; h, F. R. Goss, J. Chem. Soc. 1940, 
755. 
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From a comparison of the deviation factors obtained from the Debye relation and listed 
in the last column of table 4, it is evident that both the Onsager relation and the present modi- 
fication lead to a marked improvement with regard to the apparent variation in dipole moment 
from solvent to solvent. This is most evident in the four prolate cases for which values of 
(j Debye decrease from 20 to 35 percent in going from media of lowest to highest dielectric con- 
stant, whereas there is little or no trend in the values of G e . 



Table 4. 



Hedestrand coefficients and G factors for dilute solutions of several polar solutes in various nonpolar 

solvents 



Solute 



C 2 H 5 Br *.. 
(20° C). 



w-CWh- 
C 6 Hi2-.- 
CC1 4 _ — 
C 6 H 6 -_-_ 

cs 2 



C 6 H 5 Bra S 6 Jj 12 --- 

(20° o gSt::: 

cs 2 



Solvent 



CH3CN b_. 

(25° C)-. 



C 6 H 5 CN ■„. 
(20° C)~. 



(CH 3 ) 3 CC1 *>_. 
(25° C) 



CHCI3 <*- 
(20° C). 



w-C 6 Hh-- 
Dioxane.. 

CCI4 

C 6 H 6 

cs 2 

Tetralin. 

n-CtELu.. 

C 6 H,2...- 

CCI4 

CeH 6 

CS 2 



n-C«Hi4-- 

Dioxane- 

CCI4 

C 6 H 6 

CS 2 

Tetralin.. 



»-C»Hi4_ 

CioH, 8 ... 
CC1 4 _ — 
C 6 H 6 — . 

CS 2 



1.90 
2.02 
2.24 

2.28 
2.64 

1.90 
2.02 
2.24 
2.28 
2.64 

1.88 
2.20 
2.23 
2.27 
2.63 
2.75 

1.90 
2.02 
2.24 

2.28 
2.64 

1.88 
1.91 
2.20 
2.23 
2.27 
2.63 
2.75 

1.90 
2.20 
2.24 
2.28 
2.64 



dtld/2 



3.18 
3.85 
4.62 
5.03 
6.68 

2.51 
2.90 
3.26 
3.51 
4.34 

8.75 
17.5 
14.5 
16.5 
23.9 
10.3 

14. 16 
17.25 
21.19 
22.97 
33.51 

3.70 
3.33 
6.51 
5.64 
5.93 
8.29 
3.65 

1.30 
1.05 
1.64 
1.90 
2.18 



dd/df 2 



0.445 
.453 

-.088 
.488 
.225 

.692 
.67 
-.098 
.722 
.372 

.0112 
-.1428 
-.4308 
-.0466 
-.4711 
-.0640 



.188 

-.603 

.135 



.143 
.110 
-.227 
-.837 
-.038 
-.792 
-.084 

.507 

.307 

-.101 

.549 



Gon 



0.95 
.93 
.95 
.94 
.90 



.81 
.80 
.79 



.81 
.79 

.77 
.75 



.80 
.83 
.82 
.75 

1.04 
1.02 
1.15 
1.12 
1.08 
1.03 
1.04 

1.40 
1.41 
1.41 
1.47 
1.35 



1.00 
.98 
1.02 
1.01 
.97 



1.00 
1.00 
.93 

.77 
.95 
.89 
.87 
.86 
.84 

.94 
.93 



.91 

.95 
.94 
1.03 
1.01 
.97 
.90 
.90 

1.21 
1.18 
1.15 
1.21 
1.07 



ODebye 



0.97 

.92 
.94 
.90 
.76 

.90 
.90 
.85 
.81 
.67 

.71 
.80 
.75 

.77 
.65 
.60 



.86 
.82 
.71 

1.04 
1.04 
1.08 
1.05 
1.01 
.90 
.87 

1.26 
1.40 
1.41 
1.43 
1.24 



» E. G. Cowley and J. R. Partington, J. Chem. Soc. 1937, 130. 

b J. W. Smith and L. B. Whitten, Trans. Faraday Soc. 47, 1304 (1951). 

<= E. G. Cowley and J. R. Partington, J. Chem. Soc. 1936, 1184. 

d T. G. Scholte, Rec. trav. chim. 70, 50 (1951). 

Values of G 0n . and G e for the polar constituent in a number of mixtures with CC1 4 and 
other nonpolar solvents are plotted as a function of concentration in figure 3. Slight differ- 
ences between the values of G shown for the pure liquid and those listed in tables 2 and 3 are 
due to slightly different values of the dielectric constant used in the two cases. In general, 
the values of G 0n . show a definite dependence on concentration, being closer to unity for the 
dilute solutions than for the pure liquid. On the other hand, the values of G e are substantially 
independent of concentration in most cases. For the 10 systems containing an alkyl or aryl 
halide in CC1 4 , only for CH 3 I does G e differs from unity by more than 4 percent over the entire 
concentration range. 

5. Conclusion 

The preceding analyses show that a decided improvement in the original Onsager treat- 
ment is obtained when the molecular cavity is altered to accommodate nonspherical molecular 
shapes, provided the volume of the cavity is kept independent of temperature. The behavior 
of pure polar liquids over an extended range of temperature and of dilute solutions and mix- 
tures with a nonpolar solvent is, for the most part, satisfactorily correlated with a single set 
of molecular quantities involving size, shape, polarizability, and dipole moment. At the 
present time, however, limited knowledge concerning the appropriate values of polarizability, 
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size, and shape would appear to preclude the determination of dipole moments from data on 
liquids and mixtures with an accuracy comparable with direct measurements on the vapor. 
In view of the rather arbitrary assignment of cavity volume, inadequate allowance of atomic 
polarizability, neglect of optical anisotropy, and other factors that could lead to an overly 
complicated molecular model, the exact values derived for the eccentricities are not of partic- 
ular significance and should be regarded only as useful molecular parameters. 

Liquids subject to the influence of hydrogen bonding have not been considered. They 
usually depart drastically from Onsager's theory and obviously require more specialized treat- 
ment as provided by the theory of Kirkwood. This need is also supported by a variety of 
independent evidence. On the other hand, much of the evidence supporting the concept of 
close-range structure or association in the systems under consideration results from a rather 
restricted interpretation of the dielectric behavior. Usually inherent is the assumption that, 
in the absence of association, the dipole moment enhanced by its reaction field would be the 
same as obtained from Onsager's model. With the modifications proposed here, it does not 
appear necessary, in general, to take into account specific molecular interactions and the 
approximation of a continuum for the molecular environment appears adequate. 
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7. Appendix 

It is required to find the reaction field, R^, of a point polarizable dipole located at the 
center of a cavity. The cavity is imbedded in an infinite homogeneous and continuous medium 
having a macroscopic dielectric constant of e . 

The potential inside the cavity can be written 

^>=?f> +*/*>, (1) 

where <pf 0) reduces to the potential of a dipole for points in the neighborhood of the origin and 
(p( R) remains finite at all points within the cavity. Then 

R,= -V<p} R \ (2) 

For the cases that follow it is convenient to consider the corresponding problem for a 

charge, Q, situated at a point on the axis of symmetry and then to calculate ^, (m) and i? M by the 
relations 

£,= -V^ (B) . (4) 
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7.1. Sphere 

If a charge Q is located at the point (x ,0,0) then the potentials inside and outside the sphere 
are 

<Pi=Qit-%Pn(cos d) + ±B n r n P n (cos 6) (5) 

W-?S-C7,^- l P.(ooBfl5, (6) 

where (r,$,(p) are the spherical polar coordinates of a point P, and the P n 's are the Legendre 
polynomials. 

The usual boundary conditions over the surface of the sphere of radius r suffice to deter- 
mine the coefficients B n and C n . Only the B n are of interest. Calculation gives 



so that 



^ (K) -lim » r /^A_-M 2(6o-l) ,. 

Q^S* \b^J- ^ (260+1) X (8) 

7.2 Prolate Spheroid 
It is convenient to use prolate spheroidal coordinates (£,*?,£>) for which 

X=C%r] 



\^vSl > 



2/=cV($ 2 -l)(l-^ 2 )cos^ -l<n^l> (10) 

s=cVte 2 -l)(l-^ 2 )sin <? ^ ^ ^ 2tt ) 

If Z> is the distance between a point P' on the axis of symmetry (between the foci) and the 
point P, Hobson's expansion for D~ l reduces to 

1 1 00 

^=7S(2»+l)P»(u)<2»tt)/ > ,(0, (11) 

in which P nj Q n are the Legendre functions of the first and second kind. 
The required potentials can be written in the form 

^=?S(2Hi)^«ft-(ap^ / )+s«AW^(a (12) 

C k=0 k=0 

W)=-SftP*(u)<2*({). (13) 

60fc = 



The boundary conditions on the spheroid 



€l dt €j da 

give 



!. 



a k -- („ 1) (2*+ 1) Qt (MF t <&- l &toPM* U6) 



in which the prime on the function denotes differentiation. 
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The first term in the potential of the reaction field of a dipole at the origin reduces to 
and to the same approximation the reaction field of the dipole is 



(17) 



where the eccentricity e = ^~ l —^\ — (bla) 2 and 6/a is the axial radio of the spheroid. 
The reaction field coefficient previously defined becomes 



R? 2 c 



l» 



The indefinite 1 factor rVr' can be eliminated through the natural assumption thai the volumes 
of the Onsager sphere and the spheroid are equal. This condition gives 

(19) 

IT 6 

so that lii vail v 



3 26Q+1 1-g 2 
~2 € e 3 



^.fm£_L yw . (20) 



A family of curves representing k R is shown in figure 1 . 
The modified Onsager equation can be writ leu 



\t(~Z) at {e- £!kr {m.E^)du 



(21) 



in which E , m, iV, V are the applied field, molecular moment, Avagadro's number, and molar 
volume, respectively. An extension of Onsager's treatment of the sphere gives, to a high 
approximation 

£=-&Eiocai) (22) 

iiocai^ T - iL + ^L_ + ^A_. (23) 

1— «tf« 1— °Wll i- a ±^± 

in which i£ c is the cavity field, a. is the polarizability, and g=R/fjL. The subscripts || and J_ 
denote the directions parallel and perpendicular to the dipole axis. Hence, to a good approxi- 
mation 

(m>E\ r ( 1 pi , 2M /1-^/nM, m 2 1 ) , , N 

in which 



X= 



H' 






M= -zs± A v =J o 



ds 



Ze 



(s+6 2 ) 2 (s+a 2 ) !/:; 
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If the anisotropy of polarizability and shape are such that a||<7ii=<*xf7j. then 



(/l ,2M \ 

1 _4,iv K3 a "+3X a V g 
3F ( (l-a„^i) "*"& 



-iU ( 25 > 



3tr(l-« llS r,|) 
and if the polarizability is isotropic then 

The determination of #j_ follows the method used to find g\\. In the present case, the com- 
plete expansion of D~ l is required. Hobson's expansion, for P' on the minor axis (2), is 

+2S S (2»+l)(-ir \%^Ypz(v)Q:(&P7(G)PZ(I;') cos ro(*-*/2), (27) 
in which P%(x), Q%(x) are the Associated Legendre functions. The potentials are 

W Q) = <2i+S«.^»W^»(Q+2S a nm P:( v )P^(Q cos m(^-T/2) (28) 

^'=^^,w(!,(9+si:L^we:(a cos mb-*®. (29) 

€071=0 W = lTO=l 

The boundary conditions over the spheroid £=£0 yield two sets of equations to determine 
(a w , j8 n ) and (a wm , J8 WTO ). Only (a n , a ww ) are of interest and algebraic solution gives 

«,-- 2(2»+l)(e -l) Q Bfe) p :&) _ eo(2: ( &) P nfe )- (30) 

= - < 2 2 r2«4-iv ivf ^-^M 'fx n &(60<3r'(60Pir'(o)PrftO f oiN 

« MM - 7 2(2«+l)(-l) j^^J («„-!) (j^^^^- (3D 

The only terms contributing to the potential of the reaction field for a dipole at the origin 
are those containing the a nm . Hence the leading term is 



so that 



, «» _ 3 /« , ,s Q\ (60 Q\' fe)Pi (v)P\ (Q sin v 



w-l£( n Q) (60 ffi (60 ,,« 

u " ~ 2 c 3 Uo 1; <K &) PI' &0- «o<2!' (60 -Pi (60* ^ 



After algebraic reduction, the required factor becomes 

r 1+ _2_-|fk=i)_n 

L _i L ViJL «o xj 

L 60 T z-iJ 



*-»"■ ^ 
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The cavity field coefficient previously defined is 

k c =EyiE?= 2 -^L=$k R (35) 

A family of curves representing k c is given in figure 2. 

Calculation shows that for practical ranges of e the quantity 1 — a<7n/l — ag± is only slightly 
dependent upon e (and therefore upon the temperature) and except for high eccentricities and 
exceptional values of a/r 3 differs from unity by less than 10 percent. Moreover, the quantity 
(\"\~zMIL)k e differs insignificantly from unity except for very large eccentricities, so that 

°*-H(I + 1t[ 1+ i=S>"" m 

can be taken as virtually constant. 

The final equation used in the numerical computations reduces to 



e- 



47riVf a j? k c ) 3e f 

' i_ 3 V\l-ag^3kT (\-ag^2e+l KM) 



7.3. Oblate Spheroid 
The appropriate coordinates are the oblate spheroidal coordinates for which 
x=c£rj sin (p £^1 



l 



2/=cV(S 2 -l)(l->? 2 ) -l^v^y ) (39) 

z=ci,t\ cos <p 0^(p^2x '. 

Hobson's expansion for D~\ in which the point P' is on the y-axis, reduces to 

±=±T,(2n+i)P n (DQ n (iy)P n (iy'), (40) 

LJ C n = o 



in which x = ^\—7f and y = ^£ 2 —l. 
The required potentials are 

<p}v=Q±(2k+l)P k (x)QSy)PS¥)+i^ (41) 

C k=0 k=0 

ti Q) =-1bMSy)Pk(x). (42) 

*0 k=0 

The boundary conditions over the spheroid £=£ determine the (a k ,fi k ). Only terms containing 
a k contribute to the potential <p} R) . Algebraic solution gives 

^=Q (€o _ !) (2k+ !) ^^ (43) 

so that the first term in the potential of the reaction field of the dipole at the origin is 

riP"*—? ('o- i)Q»<»)Pi {ly)Pi ®. (44) 
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The reaction field of the dipole to the same approximation is 



(..- 1) [~ 1+ V L r-' t,1 '"'vF?][^ 1 -«'-"°- 1 v T^] 



C 3 €o 



in which e=£ _l - 

The reaction field coefficient, on the assumption that the volume of the spheroid and 
equivalent Onsager sphere are equal, reduces to 



l Vl-e 2 |_ « Vl-e 2 JL Vt-e 2 J 

L «o \ e vi-eVJ 






The calculation of the quantity (m-£ , )ave parallels that carried out for the prolate case. In 
view of the results obtained for this model, it will be assumed that to a good approximation 

(jjoWJ a (kgjW __1 Ijfcf (47 ) 

in which 

I _1 2ao Uo 1;A * ^ Jo (s+a 2 ) 2 (s+& 2 )>' 2 

M ■ 2 6o leo ' v " Jo (s+a 2 )(s+P) 3/2 

The cavity field coefficient is 



k c =0k B =^°±± M, (48) 



where 



r2?L~ 1+ V- tan vT^JL (1 "" 6V ^~ tan vi=^i 



(49) 



A family of curves representing k c is given in figure 2. 
The complete equation now becomes 

iV( a(|+|£/M) m 2 1 U 3e 

Calculation shows that k c (\-\-%LjM) is essentially unity for practical ranges of « and e, so that 
the working equation simplifies to 

, , , Nj a_ m 2 h ) 36 

In the case of the prolate spheroid, for which the calculations have been made, it should be 
noted that the inclusion of the correction factor 1 — ag.Jl — ag ± in the optical term of the polariza- 
tion in part nullifies the effect of the correction for shape and increases this term to a value more 
nearly that of the Onsager approximation for a sphere. Calculations for a number of specific 
examples indicates that if this correction is made, the values of G e listed in the tables will be 
reduced by about 1 percent or less. 

Washington, February 19, 1954. 
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